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A simple multiscale approach to the diffusion-driven adsorption from a solution to a solid surface is 
presented. The model combines two important features of the adsorption process: (i) the kinetics of 
the chemical reaction between adsorbing molecules and the surface; and (ii) geometrical constraints 
on the surface made by molecules which are already adsorbed. The process (i) is modelled in a 
diffusion-driven context, i.e. the conditional probability of adsorbing a molecule provided that the 
molecule hits the surface is related to the macroscopic surface reaction rate. The geometrical con- 
straint (ii) is modelled using random sequential adsorption (RSA), which is the sequential addition 
of molecules at random positions on a surface; one attempt to attach a molecule is made per one 
RSA simulation time step. By coupling RSA with the diffusion of molecules in the solution above the 
surface the RSA simulation time step is related to the real physical time. The method is illustrated 
on a model of chemisorption of reactive polymers to a virus surface. 

PACS numbers: 68.43.-h, 87.15.Rn 



Random sequential adsorption (RSA) is a classical 
model of irreversible adsorption (e.g. chemisorption) 
Given a sequence of times k — 1,2,3,..., an attempt 
is made to attach one object (e.g. a molecule) to the sur- 
face at each time point t = tk- If the attempt is successful 
(i.e. if there is enough space on the surface to place the 
molecule), the object is irreversibly adsorbed. It cannot 
further move or leave the structure and it covers part of 
the surface, preventing other objects from adsorbing in 
its neighbourhood (e.g. by steric shielding in the molec- 
ular context). 

In the simplest form, RSA processes are formulated as 
attempting to place one object per RSA time step, ex- 
pressing the simulation time in units equal to the number 
of RSA time steps k rather than in real physical time tk . 
Such an approach is useful to compute the maximal (jam- 
ming) coverage of the surface. To apply RSA models to 
dynamical problems, it is necessary to relate the time of 
the RSA simulation k and the real time tk. This is a goal 
of this paper. We consider that the adsorbing objects are 
molecules which can covalently attach to the binding sites 
on the surface. We couple the RSA model with processes 
in the solution above the surface to study the irreversible 
adsorption of molecules in real time. The time (tk — ifc-i) 
between the subsequent attempts to place a molecule is 
in general a non-constant function of k which depends 
on the kinetics of the chemical reaction between the ad- 
sorbing molecules and the surface, and on the stochas- 
tic reaction-diffusion processes in the solution above the 
surface. We illustrate our method with an example of 
the chemisorption of reactive polymers to a virus surface 
0, Sj- Finally, we show that the stochastic simulation 
in the solution can be substituted by a suitable deter- 
ministic partial differential equation which decreases the 
computational intensity of the algorithm. We show that 
it is possible to get the values of tk without doing exten- 
sive additional stochastic simulations. 

We consider a three-dimensional cuboid domain L x x 



L x x L z in which molecules diffuse (see Fig.[T|). The 
side z = of area L x x L x is assumed to be adsorb- 
ing, i.e. containing binding sites to which molecules can 
covalently attach. Our goal is to couple RSA on the 
side z — with stochastic reaction-diffusion processes 
in the solution above the adsorbing surface. Since those 
molecules which are far from the surface will have lit- 
tle influence on the adsorbtion process, it is a waste of 
resources to compute their trajectories. We will there- 
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FIG. 1: Three-dimensional cuboid domain. 

fore carefully truncate our computational domain to that 
which is effectively influenced by the reactive boundary 
at z = 0, which we denote by z < L(t). Note that L(t) 
is not fixed but a function of time — the formula for it 
will be derived later. Suppose that there are N(t) dif- 
fusing molecules in the cuboid domain L x x L x x L(t). 
Let us denote the ^-coordinate of the i-th molecule by 
Zi(t), treating molecules as points in the solution in what 
follows. Choosing a time step At, we compute Zi(t + At) 
from Zi(t), i = l,..., N(t), by 



Zi(t + At) = Zi (t) + y/2Di At 



(1) 



where £j is a normally distributed random variable with 
zero mean and unit variance and Di is the diffusion con- 
stant of the i-th molecule. In principle, we should model 
the behaviour of molecules as three dimensional random 
walks in the cuboid domain L x x L x x L(t), i.e. there 
should be equations analogous to {T]) for the x and y 
coordinates too. However, we can often assume that 
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L(t) 3> L. x in applications. Choosing the time step At 
large enough that a molecule travels over distances com- 
parable to L x during one time step, we can assume that 
the molecules are effectively well-mixed in the x and y 
directions on this time scale. Consequently, the x and 
y coordinates of molecules do not have to be simulated. 
If the original adsorbing surface is large, one often mod- 
els by RSA only a representative part of it, i.e. a square 
L x x L x which contains a relatively large number of bind- 
ing sites, but still satisfies L x -c L(t). The diffusion of 
molecules (JTJ) is coupled with other processes in the so- 
lution and on the surface as follows. 
Chemical reactions in the solution: Our illustrative 
example is the polymer coating of viruses [1, 0] • In this 
case, the polymer molecules have reactive groups which 
can covalently bind to the surface. The reactive groups 
also hydrolyse in solution. Assuming that there is one 
reactive group per polymer molecule (such a polymer is 
called semitelechelic), we have effectively one chemical 
reaction in the solution - a removal of the reactive poly- 
mers from the solution with rate A Assuming that 
XAt <C 1, the stochastic modelling of the process in the 
solution is straightforward. At each time step, the i-th 
molecule moves according to ([T]). We then generate a 
random number r; uniformly distributed on the interval 
[0,1]. If ri < XAt, we remove the molecule from the 
system. More complicated reaction mechanisms in the 
solution can be treated using stochastic simulation algo- 
rithms which have been proposed for reaction-diffusion 
processes in the literature @, @, 0]. In our case, we treat 
diffusion using the discretized version of Smoluchowski 
equation {]]). Consequently, we can follow Andrews and 
Bray [5j to introduce higher-order reactions to the sys- 
tem. 

Adsorption to the surface: The surface L x x L x at 
z = is assumed to be adsorbing. We use a simple ver- 
sion of the RSA model from [3] which postulates that the 
binding sites on the surface lie on a rectangular lattice. 
Binding a polymer to a lattice site prevents the bind- 
ing of another polymer to the neighbouring lattice sites 
through steric shielding, i.e. we consider RSA with the 
nearest neighbour exclusion as a toy model of adsorp- 
tion [![. Such a RSA model can be simulated on its own 
as shown in In this paper, we simulate it together 
with the z-variables of molecules in the solution |T]) to 
get the RSA evolution in real physical time. Whenever 
a molecule hits the boundary z = 0, it is adsorbed with 
some probability, and reflected otherwise. This partially 
adsorbing boundary condition is implemented in the RSA 
context as follows: 

(a) If Zi(t+ At) computed by |T]) is negative then, with 
probability P^fAt, we attempt one step of the RSA al- 
gorithm with the i-th molecule. If the «-th molecule is 
adsorbed, we remove it from the solution. Otherwise, 
we put Zi(t + At) = -Zi(t) - y/2D { At 



(b) If Zi(t + At) computed by Jl} is positive then, 
with probability exp[-Xi(t)xi(t + At)/(DAt)]PVAt, 
we attempt one step of the RSA algorithm with the 
z-th molecule. If the i-th molecule is adsorbed, we 
remove it from the solution. 

Here, P is a positive constant which can be related to 
the rate constant of the chemical reaction between the 
binding sites on the virus surface and the reactive groups 
on the polymer |8[ . This relation depends on the stochas- 
tic model of diffusion and for equation (TTJ) is given later 
- see formula ([5]). Conditions (a)-(b) state that only 
the fraction P^/At of molecules which hit the boundary 
have a chance to create a chemical bond (provided that 
there is no steric shielding). Obviously, if Zi(t + At) com- 
puted by (dJ is negative, a molecule has hit the boundary. 
This case is incorporated in (a). However, Andrews and 
Bray Q argue that there is a chance that a molecule 
hit the boundary during the finite time step At even if 
Zi(t + At) computed by |T]) is positive; that is, during the 
time interval [t, t + At] the molecule might have crossed 
to Zj negative and then crossed back to Zi positive again. 
They found that the probability that the molecule hit the 
boundary z = at least once during the time step At is 
exp[-Zi(t)zi(t + At)/(DAt)] for Zl (t) > 0, Zi(t + At) > 0. 
This formula is used in (b). 

Numerical results: It is important to note that the 
boundary conditions (a)-(b) can be used for any RSA 
algorithm and for any set of reactions in the solution. 
To apply it to the virus coating problem, we have to 
specify some details of the model. First of all, it can be 
estimated that the average distance between the bind- 
ing sites is about 1 nm 2]. We choose L x = 100 nm. 
Therefore, there are about 10, 000 binding sites on the 
adsorbing side 2 = 0. We use RSA on a 100 x 100 lattice 
with the nearest neighbour exclusion, which is a special 
case of the model from Q . We consider a monodisperse 
solution of semitelechelic 50 kDa polymers, i.e. Di = D 
where D = 5 x 10~ 5 mm 2 s _1 9]. The rate of hydrolysis 
of the reactive groups on polymers can be estimated from 
data in 3 as A = 1.3 X 10 -4 s _1 . We choose P = 1 s" 1 / 2 . 
Since we simulate the behaviour of polymer molecules 
in solution only along the z-direction, we express the 
concentration of polymer c(z,t) in numbers of polymer 
molecules per volume L x x L x x [1 mm] where L x = 10~ 4 
mm is fixed. A typical experiment starts with a uniform 
concentration of reactive polymers. Considering that the 
initial concentration of 50 kDa polymer is 0.1 g/1, we ob- 
tain the initial condition Co = 1.2 x 10 4 molecules per 
mm of the height above the surface (let us note that the 
units of the "one-dimensional" concentration c(z,i) are 
molecules/mm because L x is considered fixed). Next, we 
have to specify L(t) (see Fig.HJ, i-e- we want to find the 
region of the space which is effectively influenced by the 
boundary condition at z = 0. To that end, we note that 
the concentration c(z,t) satisfies the partial differential 
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equation 



dc 
Iff 



D 



dz 2 



Ac. 



(2) 



Now any partially reacting boundary will have less im- 
pact on the spatial profile of c(z,t) than perfect adsorp- 
tion at z — 0. Thus we may find an upper bound for the 
region of influence of the boundary by solving ((2|) subject 
to 

c(0,t) = 0, lim c(z,t) = c exp[— At], (3) 

2 — >00 

for t G [0, oo), and the initial condition c(z,0) = cq, for 
z e [0,oo). The solution of ©-© is 



c(z, t) — cq exp[— At] erf 



2VDi 



(4) 



where erf(-) denotes the error function. Defining w 
erF 1 (0.99) = 1.821 we set 

L(t) = 2ujy/Dt. 



(5) 



Then c(L(t), t) = 0.99 c exp[— At], so that the concentra- 
tion of the reactive polymer at point z = L(i) at time t 
is equal to 99 % of the polymer concentration at points 
which are "infinitely far" from the adsorbing boundary. 
In particular, we can assume that the adsorbing bound- 
ary effectively influences the polymer concentration only 
at heights z S [0, L(t)] above the boundary and we can 
approximate c(z,t) ~ coexp[— Xt] for z > L{t). Formula 
(O specifies the computational domain as L x x L x x Lit) 
at each time (see Fig.QJ. 

The results of the stochastic simulation of the solution 
above the surface are shown in Fig. [2] as grey histograms. 
To simulate the behaviour of N(t) reactive polymers, we 
consider only their z-positions. We use At = 10 -2 s and 



we update the z-positions of molecules during one time 
step according to ([J). At each time step, we also gen- 
erate a uniformly distributed random number and we 
remove the z-th molecule from the system if < A At. 
We work in the one-dimensional domain [0, L(t)] where 
L(t) is given by (JSJ). The RSA boundary condition at 
z = is implemented using (a)-(b) described above. 
The right boundary increases during one time step by 
AL(t) = L(t+At) — L(t). During each time step, we have 
to put on average m(t) = Cq exp[— At] AL(t) molecules 
into the interval [L(i),L(t + At)]. This is done as fol- 
lows. We put L m (^)J molecules at random positions in 
the interval [L(t),L(t + At)], where |_ - J denotes the in- 
teger part. Moreover, we generate random number tm 
uniformly distributed in [0, 1] and we add one molecule 
at a random position in the interval [L(t),L(t + At)] if 
fAt < m(t) — \m(t)] . This will ensure that we put on 
average m(t) molecules to the interval [L(t),L(t + At)] 
during one time step. 

Introducing the moving boundary decreases the com- 
putational intensity of the model. Initially we simulate a 
relatively small region with a high concentration of reac- 
tive polymers. The simulated region increases with time 
but the concentration of reactive molecules decreases 
with rate A. Using |^|). it can be computed that the maxi- 
mal number of simulated polymers in solution is achieved 
at time t m — (2A) _1 = 64 min (and is about 8 x 10 3 
molecules for our parameter values). 

The number of polymers adsorbed to the RSA surface 
at z = as a function of real physical time is shown 
in Fig.[3l Since the polymer solution is assumed to be 
monodisperse, we can run the RSA algorithm first and 
record the times ki , It2, k^, ... (expressed in numbers of 
the RSA time steps) of successful attempts to place the 
polymer on the RSA lattice. Then the stochastic simula- 
tion of the reaction-diffusion processes in the solution can 
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FIG. 3: Number of polymer molecules adsorbed to the RSA 
lattice as a function of the real physical time. 



use ki, k%, &3, . . . as its input. We will shortly consider 
another approach to the problem, replacing the stochas- 
tic simulation of the solution by the continuum limit ([2|) 
with a suitable Robin boundary condition. To enable 
a direct comparison of the two approaches, we use the 
same sequence k\ , £>2 > &3 , ■ ■ ■ in ten realizations of the 
full stochastic model of adsorption; the results are shown 
as grey solid lines in Fig.[3J 

RSA-PDE approach: Moving the right boundary L(t) 
is one way to decrease the computational intensity of the 
problem. Another possibility is to use the deterministic 
equation |J3J) together with a Robin boundary condition 



dc 

dz 



(0,i) = a RSA (t) 



2P 



c(0,t) 



(6) 



at the adsorbing boundary z — 0. Here, the fraction 
2P/tJDit corresponds to the rate of the chemical reac- 
tion between the adsorbing boundary and the diffusing 
molecules in the solution - see [H for the derivation of this 
formula and further discussion. Factor aRSA{t) £ {0,1} 
provides the coupling between the RSA model and ([2]). 
To find the value of ctRSAit), we estimate the number of 
attempts to place the polymer on the RSA lattice by 



K(t) 



2PVD 



c(0,t)dt 



where |_ - J denotes the integer part |8|. We start with 
&RSA{t) — 1 and we solve © and © numerically. When- 
ever n(t) increases by 1, we attempt one step of the RSA. 
If the attempt is successful, we put ansA(t) = 1. If the 
attempt to place the molecule is not successful, we put 
&RSA{t) = 0. Thus ctRSA{t) has only two values, and 1, 
and changes at computed time points depending on the 
output of the RSA simulation. We call this procedure 



the RSA-PDE approach. It also leads to the sequence of 
real physical times t] il , tk 2 , ifc 3 . . . , of successful attempts 
to place the polymer on the RSA lattice. The numerical 
solution of equation {2J with the Robin boundary condi- 
tion © at z = is presented in Fig.f2]as the solid line 
for comparison. We also plot the number of adsorbed 
polymers as a function of the real time as the dashed line 
in Fig. [3] To enable the direct comparison, we run the 
RSA algorithm first and we record the times of successful 
attempts to place the polymer on the lattice. We obtain 
the sequence fci, &2, &3, ... of times expressed in number 
of RSA time steps. This sequence is used in both the 
stochastic model (10 realizations plotted in Fig. [3] as grey 
solid lines) and the RSA-PDE approach (dashed line in 
Fig. [3]). The comparison of the results obtained by the 
full stochastic model and by the RSA-PDE model is ex- 
cellent. 

Conclusion: We have presented a method to perform 
RSA simulation in real physical time. The key part of 
the method is the boundary condition (a)-(b) which can 
be coupled with any reaction-diffusion model in the solu- 
tion and any RSA algorithm. We illustrated this fact on 
a simple model of the polymer coating of viruses. More- 
over, we showed that the RSA algorithm can be coupled 
with @ using the Robin boundary condition ^ to get 
comparable results. The Robin boundary condition (J6]) 
is also not restricted to our illustrative example. It can 
be used for the coupling of any RSA model with the PDE 
model of the reaction-diffusion processes in the solution 
above the adsorbing surface. 
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